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Abstract: The six-fold rotational symmetry of photonic crystal fibers has 
important manifestations in the radiated fields in terms of i) a focusing 
phenomena at a finite distance from the end-facet and ii) the formation 
of low-intensity satellite peaks in the asymptotic far field. For our study, 
we employ a surface equivalence principle which allows us to rigorously 
calculate radiated fields starting from fully-vectorial simulations of the near 
field. Our simulations show that the focusing is maximal at a characteristic 
distance from the end- facet. For large-mode area fibers the typical distance 
is of the order 10 x A with A being the pitch of the triangular air-hole lattice 
of the photonic crystal fiber. 
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Fig. 1. Schematics of the fiber geometry and the imaginary closed surface S = Si -\- S2, 
where is a circle parallel and close to the end-facet and S2 is a semi- sphere concentric 
withal. 

The photonic crystal fiber [1,2] (PCF) offers a unique opportunity for realizing broadband 
large-mode area (LMA) single-mode (SM) fibers [3]. In high-power laser applications the LMA 
property is desirable to keep the power density below the damage threshold of silica as well as 
to avoid non-linear phenomena and the SM property is essential for the beam quality and the 
ability to generate close to diffraction-limited beams. Recent studies of LMA-SM-PCF lasers 
have suggested that the PCF technology may indeed offer good beam quality [4] as quantified 
by an M-squared value [5] close to unity (M^ = 1 for a Gaussian beam). 

For the fundamental mode in PCFs the near and far fields have intensity distributions which 
are overall close to Gaussian, but with a six-fold rotational symmetry reflecting the symmetry 
of the air-hole lattice. For LMA-PCFs there are no significant deviations from a Gaussian dis- 
tribution for intensities above the l/e'^ level [6] and for many practical applications the mode 
may be considered Gaussian. However, as studied recently the derivation from a Gaussian dis- 
tribution has some important consequences [7] . In the near-field the intensity is suppressed at 
the locations of the air-holes and in the asymptotic far-field these suppressions manifest them- 
selves as low-intensity satellite peaks. Furthermore, in the near field to far-field transition the 
intensity becomes strongly focused at a finite distance from the end-facet as can be understood 
by a simple Gaussian decomposition of the near field. 

While the simple decomposition approach has proved very useful for the qualitative under- 
standing of the near to far-field transition the laser applications call for a more quantitative study 
of the radiated field. Calculating the diffraction pattern is generally a complicated task, which 
involves the full solution to the elastic scattering problem at the end-facet. Here, we rigorously 
calculate radiated fields from the end-facet utilizing a surface equivalence principle [8]. 

The calculations of the modal reflectivity from the fiber end-facet is performed with a three- 
dimensional frequency domain modal method relying on Fourier expansion technique [9, 10]. 
In brief, the Fourier expansion method relies on an analytical integration of Maxwells equations 
along one direction (the fiber axis) and on a supercell approach in the other two transversal di- 
rections. We emphasize that for the radiated field at the surface Si we have treated the elastic 
scattering at the end-facet fully, i.e. the incident fundamental mode is partly reflected (no mode 
mixing for symmetry reasons) and partly radiated/transmitted. Through the surface equivalence 
principle, any wave-front can be considered as a source of secondary waves which add to pro- 
duce distant wave-fronts according to Huygens principle. 

In order to apply the surface equivalence principle, we surround the end-facet radiating fiber 
with an imaginary closed surface S = Si -\- S2, where Si is a circle parallel and close to the 
end-facet and 5*2 is a semi-sphere concentric with Si, see Fig. 1. Both of these surfaces are of 
infinite extent. By considering proper surface currents on S which are equivalent in the sense 



that they yield the same field outside S, we can formally erase all the existing elements inside 
S and fill it with a perfect electric conducting (PEC) body. The equivalent surface electrical 
and magnetic currents on S are Je(r) = n x H and Me(r) = — n x E, respectively, where n is an 
outward normal to the surface S [8]. Of course, these currents radiate in the presence of the PEC 
and not in unbounded free- space. Because of the infinite radius of ^2, electromagnetic fields on 
this surface take zero value. On Si electromagnetic fields are determined from near-fields. Since 
the surface is an infinitely extended flat PEC body we can utilize image theory to replace 
the conductor by image currents — Je and M^. Ultimately, computation of radiated fields of 
an open-ended PCF is simplified to computing radiated fields of the fictitious current source 
M(r) = 2Mg(r) = — 2n x E in free-space. Calculating radiated fields of the current source can 
be achieved systematically through electric 7re{r) and magnetic 7r^(r) Hertz vector potentials 
defined by: 



where CO = ck, k = In/A., and A is the free-space wavelength. From Eqs. (2a) and (2b) we 
expect an improvement in computational time by a factor of roughly two compared to a direct 
use of the surface equivalence principle [11] in which both fictitious electric and magnetic 
surface currents exist. 

We employ the described approach in a fully vectorial study of the evolution of the funda- 
mental mode of a PCF from near-field to far-field. As an example we consider the fundamental 
mode in a pure-silica PCF with a normalized air-hole diameter d/A = 0.45 and in our sim- 
ulations we utilize the scale-invariance of Maxwell's equations for a frequency-independent 
dielectric function. Recent experimental studies showed quite complicated interference phe- 
nomena in the near to far-field transition [7] and in particular a focusing behavior was observed 
at a finite distance from the end-facet. As illustrated in Fig. 2 this focusing phenomena is born 
out by our numerical simulation where the peak intensity of the radiation increases up to short 
distances from the fiber end-facet. The electric field intensity is maximal at a distance zo which 
is typically in the range A < zo ^ 10 x A. At this particular distance the intensity pattern is 
rotated by an angle of n /6 compared to the near-field, in full agreement with the experimen- 
tal observations. As the wavelength is increased relative to the pitch the diffraction by the six 
inner-most air holes increases and zo shifts closer to the end facet of the PCF. We emphasize 
that this interference phenomenon has no counterparts for e.g. step-index fibers with cylindrical 
symmetry; the focusing relies fully on the six-fold rotational symmetry of the air-hole lattice. 
It should also be emphasized that the phase-front is non-trivial (non-constant) in the transverse 
plane at zo and in that sense the focusing differs from that observed with a classical lens. 

In Fig. 3 we show the electric field amplitude at various distances from the end-facet up to 
the asymptotic far-field region. The nearly hexagonal shape of the fundamental mode of the 
fiber (Panel A) is transformed into a nearly circular shape (Panel B) which is followed by a 
hexagonal shape at z — zo (Panel C) where the orientation is rotated by an angle of ;r/6 with 
respect to the fiber mode at z = 0. This hexagonal shape expands as it propagates (Panels D and 
E) and then becomes close to circular again above the l/e amplitude level (Panel F). As the 
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Employing these potentials, electric and magnetic fields are obtained by: 



E(r) = + V^)7ie{r) - jcoiio^ x 7r^(r), 
H(r) = j(0£oV X Tieir) + {k^ + ^^)n^{r), 
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Fig. 2. Variation of the electric field intensity with distance fi-om the end-facet (z = 0) at 
the center of fiber. 



distance from the end-facet becomes larger, six satellites start emerging from the mode shape 
below the l/e amplitude level (Panel G). Finally, for larger distances satellites appear clearly 
(Panel H). 

As an alternative to the above rigorous approach, it is in fiber optics common to approximate 
the far field by a Fourier transformation of the near field [12]. Here we shall derive this expres- 
sion starting from the above formalism. In the far-field limit we approximate the |r — r'| in the 
denominators of Eqs. (la) and (lb) by r and in the exponential we correspondingly approximate 
it by r — f • r', where f = r/r is the unit radial vector of the spherical coordinate. The magnetic 
Hertz vector potential, Eq. (lb), then simplifies to 
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with 6 and being the polar and azimuthal angles, respectively, and r = rr -\- (j)^ -\- 06. The 
electric field now becomes 
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We will consider small angles of divergence so that f x M = — MyX + M^y -\- ^{0)hy which 
Eq. (5) simplifies to 
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Fig. 3. Electric field amplitudes at distances z from the end-facet varying from z = (Panel 
A) to z = 12 X A (Panel H). 

with the subscript t indicating the transverse component. The integral is nothing but a two- 
dimensional Fourier transform of the near-field corresponding to the overlap of the near field 
with a free- space plane- wave. 

Fig. 4 shows a plot of the asymptotic far-field intensity, including results of the rigorous 
surface-equivalence principle as well as the approximate far-field calculated from Eq. (6). The 
far-field pattern is computed for A /A = 0.1 at a distance of 1000 x A from the end-facet, 
which is far from the starting point of the Fraunhofer region. For an extended source of electro- 
magnetic radiation, the Fraunhofer region appears at a characteristic distance 2D^ /X from the 
source where the parameter D is defined by the maximum linear dimension of the source. For a 
PCF D ~ 2A, so Fraunhofer region sets in around 80 x A. The inset shows a contour plot of the 
far-field intensity calculated with the rigorous surface-equivalence principle. We observe six 
low-intensity satellite peaks as well as further high-order diffraction peaks, in full agreement 
with experimental observations [1, 7]. As seen, the approximate results of Eq. (6) are in quali- 
tative agreement with the full numerical solution, and it is only in the sub-1/^^ intensity regime 
that we observe quantitative deviations. This suggests that Eq. (6) could be a good starting point 
for numerical calculations of the NA and the M- squared value for LMA PCFs, which we will 
report elsewhere. 

In conclusion we have used a surface equivalence principle to rigorously calculate radiated 
fields starting from fully-vectorial simulations of the near field. Our simulations have revealed 
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Fig. 4. Asymptotic far-field intensities calculated by the rigorous approach as well as the 
approximate expression, Eq. (6), along the high- symmetry directions (j) =0 and (j) = 7i/2. 
The inset shows the corresponding full contour plot obtained with the rigorous approach. 

a focusing behavior which is maximal at a characteristic distance, of the order 10 x A from 
the end-facet. In the far-field limit we have shown how qualitative and to some degree also 
quantitative insight may be gained from simple two-dimensional Fourier transforms of the near- 
field. 
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